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ABSTRACT 

Equipartition magnetic fields can dramatically affect the polarization of radiation 
emerging from accretion disk atmospheres in active galactic nuclei. We extend our pre- 
vious work on this subject by exploring the interaction between Faraday rotation and 
absorption opacity in local, plane-parallel atmospheres with parameters appropriate 
for accretion discs. Faraday rotation in pure scattering atmospheres acts to depolarize 
the radiation field by rotating the polarization planes of photons after last scattering. 
Absorption opacity in an unmagnetized atmosphere can increase or decrease the polar- 
ization compared to the pure scattering case, depending on the thermal source function 
gradient. Combining both Faraday rotation and absorption opacity, we find the follow- 
ing results. If absorption opacity is much larger than scattering opacity throughout the 
atmosphere, then Faraday rotation generally has only a small effect on the emerging 
polarization because of the small electron column density along a photon mean free 
path. However, if the absorption opacity is not too large and it acts alone to increase 
the polarization, then the effects of Faraday rotation can be enhanced over those in 
a pure scattering atmosphere. Finally, while Faraday rotation often depolarizes the 
radiation field, it can in some cases increase the polarization when the thermal source 
function does not rise too steeply with optical depth. We confirm the correctness of 
the Silant'ev (1979) analytic calculation of the high magnetic field limit of the pure 
scattering atmosphere, which we incorrectly disputed in our previous paper. 
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1 INTRODUCTION 

Explaining the optical polarization observed in active galac- 
tic nuclei (AGN) has long been a problem for accretion disk 
models. Optically thick, pure electron scattering discs are 
expected to emit radiation which is linearly polarized up to 
11.7 percent parallel to the plane of the disk (Chandrasekhar 
1960), but this has never been observed in type 1 AGN. One 
possible reason is that the optical radiation emerging from 
the disk is Faraday depolarized by photospheric magnetic 
fields. In a previous paper (Agol & Blaes 1996, hereafter 
paper I), we have shown that such fields can drastically re- 
duce the polarization at optical wavelengths if they are near 
equipartition strength. 

Absorption opacity can also play a significant role in 
determining the polarization of the emerging radiation field. 
In a simple investigation of this effect, Laor, Netzer & Pi- 
ran (1990) showed that absorption opacity can reduce the 
overall polarization. A more careful treatment by Blaes & 
Agol (1996) showed that while this is often qualitatively 
true, absorption opacity can sometimes increase the over- 



all polarization by increasing the limb darkening (see also 
Bochkarev, Karitskaya, & Sakhibullin 1985). 

The Faraday rotation calculations of paper I assumed a 
pure electron scattering atmosphere, but absorption opacity 
might reduce the depolarization. This is because the Fara- 
day rotation of a given photon depends on the total electron 
column density that the photon traverses. The dominant 
effect of Faraday rotation occurs after last scattering (pa- 
per I) , so if the absorption opacity significantly reduces the 
electron scattering column down to unit optical depth, the 
depolarization would be smaller. On the other hand, the ab- 
sorption opacity itself may directly increase or decrease the 
polarization from an unmagnetized disk, as noted above. In 
this paper we attempt to disentangle these effects in order 
to understand how Faraday rotation and absorption opacity 
act together to determine the polarization of the radiation 
emerging from AGN accretion discs. We have discovered a 
number of subtle phenomena which are not immediately ob- 
vious from the above arguments. 

As in paper I, we use Monte Carlo calculations of the 
radiative transfer. In addition, however, we also show in sec- 
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tion 2 below how Faraday rotation by a uniform, vertical 
magnetic field can be incorporated directly into the radia- 
tive transfer equation. The emerging radiation field can then 
be calculated much faster using standard finite difference 
techniques, and we present the results of both approaches 
in simple toy atmosphere models in section 3. In section 4 
we discuss again the role of Faraday rotation in determin- 
ing the optical polarization in AGN accretion discs, and we 
summarize our conclusions in section 5. 



2 EQUATIONS AND NUMERICAL 
TECHNIQUES 

To calculate the radiation field emerging from the accretion 
disk, we treat each portion of the disk photosphere as a 
locally plane-parallel, semi-infinite atmosphere. At the op- 
tical and ultraviolet photon frequencies of interest, electron 
scattering in the magnetized plasma has negligible circular 
dichroism, so the polarization of the radiation field must be 
nearly linear (cf. paper I).^ 



2.1 Monte Carlo Calculations 

Let p be a unit vector in the plane of polarization of a given 
photon, perpendicular to its direction of propagation n, also 
a unit vector. Between scatterings, the photon polarization 
will be Faraday rotated to 



Prot = p cos X + (n X p) sin x, 



(1) 



where x = b ■ nTT5/2, tt is the Thomson scattering depth 
along the photon trajectory, and b = B/B is a unit vector 
along the magnetic field B. The photon wavelength A and 
magnetic field strength only enter through the parameter 



V5000A/ 



(2) 



where e is the electron charge. 

Paper I describes a Monte Carlo technique based on 
these equations to calculate the polarized radiative transfer 
through a magnetized, pure electron scattering atmosphere. 
We have modified this slightly to include the effects of ab- 
sorption opacity at frequency i/, assuming for simplicity 
that the ratio of absorption opacity to electron scattering 
opacity is independent of optical depth in the atmosphere. 
In other words. 



HeUT 



+ UeO-T 



(3) 



* Whitney (1991a, 1991b) has conducted Monte Carlo calcu- 
lations of the polarization of a magnetized, electron scattering 
atmosphere. Her calculations provide an interesting contrast to 
ours, because she included magnetic corrections to the scattering 
cross-section, but neglected Faraday rotation. Her results are of 
relevance to magnetic white dwarf and neutron star atmospheres. 
Our work has neglected magnetic effects on the scattering cross- 
section but has included Faraday rotation. This is much more 
relevant to optical and ultraviolet radiation emerging from AGN 
accretion discs, because the corrections to the scattering cross- 
section are negligible. 



is constant, where rie is the electron number density and ar 
is the Thomson cross section. 

We propagate each photon a vertical optical depth 
Ti, = r^o + piln(r-i) through the atmosphere, where ri is 
a random deviate between and 1, ^ is the direction co- 
sine of the photon propagation vector with respect to the 
upward vertical, and t^o is the starting optical depth. The 
photon's polarization vector is Faraday rotated according to 
equation (Q). Then, another random deviate, r2, between 
and 1 is chosen. If r2 is less than g^, the photon is scattered. 
Otherwise it is absorbed and another photon is started at 
the base of the atmosphere. This process is repeated until 
a photon escapes from the atmosphere, and it is binned as 
described in paper I. 

2.2 Feautrier Radiative Transfer 

One of the advantages of the Monte Carlo technique is that 
it is capable of handling general, complex geometries. Be- 
cause our purpose in this paper is to understand the physi- 
cal effects of Faraday rotation in the presence of absorption 
opacity, we limit consideration to locally plane-parallel at- 
mospheres with uniform, vertical magnetic field. The radia- 
tion field will then be completely axisymmetric and depend 
only on vertical depth. In this case it is straightforward to 
include the Faraday rotation directly in the full radiative 
transfer equation by just adding an extra term. This equa- 
tion can then be solved much more quickly using standard 
numerical techniques. 

The full polarized radiative transfer equation for a gen- 
eral magnetoactive plasma is already well-known (see e.g. 
Silant'ev 1979). However, because in our case it is so simple 
and illuminates the physics, we now briefly sketch a deriva- 
tion of the Faraday rotation term. 

We first project the photon polarization vector on two 
orthogonal axes which are perpendicular to the propagation 
direction n. Let the first axis be parallel to the plane of 
the atmosphere, and the corresponding polarization vector 
component be po. Let the polarization vector component 
with respect to the second axis be pgo. Then equation (|l]) 
implies that after Faraday rotation. 



Porot = Po cos X - P90 sm X 
and 

P90rot = po sin X + P90 cos X- 



(4) 



(5) 

Following Chandrasekhar (1960), define Irv and Jj^ as 
the intensities of the radiation corresponding to po and poo, 
respectively. The Stokes parameter Q,^ may them be de- 
fined as Irv — hv In a similar fashion, let Uv be the Stokes 
parameter with respect to two axes rotated by 45 degrees 
from those defined previously. The Stokes parameter Vv van- 
ishes because the radiation is linearly polarized. Expressed 
in terms of the total specific intensity 1^ = In^ + hu and 
averages over the individual polarization vectors of the cor- 
responding photons, we have = Iu{< Pq > — < Pgo >) 
and Uv = Ii, < 2poP9o >■ The degree of polarization is 



(6) 



(Note that, in contrast to paper I, we have not normalized 
the Stokes parameters by the total intensity.) 
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We now describe the intensity and polarization of the 
radiation field with the column vector 




(7) 



From equations (^) and (^), Faraday rotation transforms the 
radiation field according to 

no \ 

I.rot = cos2x -sin2x I.. (8) 
\ sin 2x cos 2x ) 

Let z be the height measured vertically upward in the at- 
mosphere. Then for an infinitesimal change in height dz, the 
corresponding Faraday rotation angle for a vertical magnetic 
field is 



(9) 



(10) 



where the upper (lower) sign is to be taken for an upward 
(downward) directed field. Expanding equation (^), we de- 
duce that the effect of Faraday rotation by a vertical mag- 
netic field can be described by 

dz 
where 

/O 0\ 

F = ±5 1 . (11) 
Vo -1 0/ 

Note that F does not change the total intensity, as expected. 

Inserting this term into the full radiative transfer equa- 
tion, we have 



91. 



r 



(12) 



Vo/ 

+ ^neCrT J dfj,'P(p, 

where r;. is the thermal emission coefficient (assumed unpo- 
larized) , 

p(m,m') = 



I [l + iP2(M)P2(M')] (1-m")J'2(m) 0^ 

(i-m')P2(m') |(i-m')(i-m") 

0, 



(13) 



and P2{fJ.) = (3/i^ - l)/2 is a second order Legendre poly- 
nomial (Chandrasekhar 1960, Loskutov & Sobolev 1979). 
Switching to the total optical depth as the dependent 
variable in the usual way, 




(14) 

where = 77. /k. is the thermal source function. 

The formal solution for equation ( |l^ can be expressed 
in terms of the total source function: 

3 = 1^ 4 ^ = S{l-q) VP(m, m')I(m')-(15) 
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(Note that = and 9q only has a contribution from 
scattering.) Then, the formal solution is given by: 



1(0, m) = 



,dt 



'10 ° \ 

costSq ±sint5g 1 Ss(f, /x)e~*'''' — ,(16) 
, =psint5g costSq / ^ 



where the sign convention is the same as for equation (g). 
The matrix represents the effect of Faraday rotation from 
the point of emission or last scattering to the top of the 
atmosphere (cf. eq. 

We have applied the Feautrier technique (e.g. Mihalas 
& Mihalas 1984, Phillips & Meszaros 1986) to solving equa- 
tion ( p^ subject to the boundary condition that there be 
no external illumination of the atmosphere at = 0. Unless 
otherwise noted, we calculate the integrals over ji with six- 
teen point Gaussian quadratures and use a logarithmically 
spaced grid in r. 



3 POLARIZATION FROM CONSTANT 
ATMOSPHERES 

In order to illuminate the physics, we now consider two ide- 
alized atmosphere problems, both with q^ independent of 
optical depth. The first case has zero thermal source func- 
tion everywhere except for a source at infinite optical depth. 
The second case has an isotropic thermal source function 
which varies linearly with optical depth. These problems 
were solved in the unmagnetized case by Loskutov & Sobolev 
(1979). Photons of different frequency are completely decou- 
pled in our radiative transfer equation for an atmosphere of 
fixed assumed structure, i.e. there is no frequency redistri- 
bution. We therefore drop the subscript v on all variables 
and parameters from now on. 



3.1 Case 1: Radiation Sources from Infinite 
Optical Depth 

This problem is a generalization of the pure electron scat- 
tering case considered by Chandrasekhar (1960), but it is 
important to note that the presence of absorption opacity 
in this case implies that the intensity of the radiation emerg- 
ing from the top of the atmosphere is formally zero unless 
there is infinite illumination from below. Nevertheless the 
radiation field is still polarized. This problem therefore rep- 
resents an idealization of an atmosphere in which photons of 
a given frequency are thermally emitted in significant quan- 
tities only at large optical depth. 

From the numerical standpoint we have performed both 
the Monte Carlo simulations and the Feautrier calculation 
using sources at sufficiently high, but finite, optical depth so 
that the polarization no longer depends on this depth. We 
apply a lower boundary condition of unpolarized, isotropic 
radiation sources. The results then depend on only two pa- 
rameters, 5 and q. 

Figure 1 shows the polarization as a function of viewing 
angle for a variety of values of q in an unmagnetized atmo- 
sphere. The lowest curve shows the pure electron scattering 
case {q — 1). Loskutov & Sobolev (1979) numerically cal- 
culated cases for q > 0.5, and found that the polarization 
increased monotonically as q decreased. Their results are 
also shown in figure 1 and are in excellent agreement with 
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Figure 1. Polarization as a function of viewing angle for an un- 
magnetized (5 = 0) atmosphere with different values of q and all 
radiation sources at infinite optical depth. From top to bottom, 
the curves represent the results of our Feautrier calculations for 
g = 0.1 to g = 1 in steps of 0.1. Points represent the numerical 
results of Loskutov & Sobolev (1979), which are consistent with 
our results. The dashed line represents our analytic solution for 
9^0. 



Figure 2. Polarization as a function of viewing angle for (a) 
q = 0.8 and (b) q = 0.2 atmospheres with various values of <5 and 
all radiation sources at infinite optical depth. From top to bottom, 
the curves represent the results of our Feautrier calculations for 
(5 = (i.e. zero magnetic field), 2, 5, 50, and 100. Square points 
represent the results of our Monte Carlo calculations for the case 
q = 0.8 and 5 = 5, and confirm the Feautrier results. 



ours. They also found that the polarization should continue 
to rise for even smaller q, and we again confirm this fact. 
Physically this is somewhat puzzling, because it suggests 
that the polarization remains finite even in the q ^ Q limit 
where there is no scattering opacity. 

We have obtained the following g ^ analytic solution 
to this problem in the unmagnetized case. For outward rays 
(0<M<1), 



and 



(1 - ^)e*/39 + 2 



/(0,1), 



Q(r,M) = (l + M)e""''''/(0,l). 
For inward rays (—1 < < 0), 

i{t,n = iy^) -'(0,1), 

and 

Q(r,M) = (1 + M)(e" - e^/^)e-*/^'/(0, 1). 



(17) 



(18) 



(19) 



(20) 



The Stokes parameter U vanishes because 5 = 0. This so- 
lution may be verified directly by substitution in equation 
(p^. In this limit the emergent intensity vanishes except in 
the upward vertical direction (/^ = 1), i.e. there is absolute 
limb darkening. The polarized fiux which is represented by Q 
vanishes for all viewing angles, consistent with the fact that 
there is no scattering opacity. The degree of polarization 
does not vanish, however, except along the vertical (/i = 1) 
because of symmetry: 



(21) 



This polarization is also plotted in figure 1 as the dashed 
line. 

For magnetized atmospheres with finite 5, equations 
(p^-(pl|) still represent the solution for the radiation field in 
the (jf — > limit. This is because the Faraday rotation term 
in equation ( [l4[ ) is proportional to q. In other words, Fara- 
day rotation depends on the electron density, and therefore 
must have negligible effect when absorption dominates over 
electron scattering. This is true even though this electron 
scattering produces a nonzero degree of polarization. 

Figure 2 shows the polarization for q = 0.8 and q = 0.2 
atmospheres with various magnetic field strengths repre- 
sented by 5. As expected, Faraday rotation depolarizes the 
radiation field. This figure should be compared to figure 
2(b) of paper I which shows the same thing for a g = 1 
(zero absorption, pure electron scattering) atmosphere. It is 
apparent that moderate absorption opacity (e.g. q = 0.8) 
enhances the depolarizing effects of the magnetic field, even 
though the electron scattering depth down to unit optical 
depth is smaller. This is true even along lines of sight which 
are perpendicular to the magnetic field (^ = 0). 

Based on our Monte Carlo results in paper I, we had 
claimed that Faraday rotation by a vertical magnetic field 
appeared to have no effect on the polarization for = 0, at 
least for the q = 1 case considered there. We rationalized 
this result on physical grounds by noting that Faraday ro- 
tation probably had its primary effect on photons after last 
scattering, and such photons would not suffer any rotation if 
they travel perpendicular to the magnetic field. However, we 
have repeated the q ~ 1 calculations with our Feautrier code 
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Figure 3. Polarization multiplied by (1+^^(5^)^''^ for various 5's. 
The solid curve represents Silant'ev's formula for large 5, while 
the dashed curves and points depict the results of our Feautricr 
code. 



and have found that the polarization at /i = drops from 
11.7 per cent at 5 = (the value from Chandrasekhar 1960) 
to 9.14 per cent as 5 — > oo. This agrees with the analytic 
calculation of Silant'ev (1979), which we wrongly disputed 
in paper I. Our Monte Carlo simulations did not have the 
resolution to see this decrease since the number of photons 
at /i = is so small. 

Figure 3 shows a comparison of our numerical results 
with the high S calculation of Silant'ev (1979). It is quite 
challenging to calculate the high S case numerically. For 
5 = 50, we reached convergence only with 20,000 logarith- 
mically spaced depth points from r = 10~^ to t = 10 for 
8 angular points. Even larger S becomes numerically pro- 
hibitive because of the large number of matrices that need 
to be stored in the Feautrier method. In any case we do find 
agreement with Silant'ev's formula for large S. (In the ap- 
pendix we show how Silant'ev's result can be derived from 
our radiative transfer equation p^ ) 

Since our Feautrier calculations are more accurate than 
the q = 1 Monte Carlo calculations from paper I, we can 
assess the accuracy of the analytic fitting formulae found in 
paper I. For the q — 1 case with 5 < 10, our fitting formulae 
are accurate to better than 12 per cent for P, 9 per cent for 
Q, and 35 per cent for U (it is least accurate when is very 
small). Silant'ev's formula can be used when 5 > 10, where 
it is accurate to better than 11 per cent for P. 

Faraday rotation generally has only a very small effect 
on the limb darkening of the total intensity of the radia- 
tion field emerging from an atmosphere, as shown in figure 
4. In the pure scattering, q = 1 case, large Faraday rota- 
tion acts to randomize a photon's polarization vector be- 
tween scatterings, and the limb darkening law approaches 
that for scattering described by Rayleigh's phase function 
(cf. the appendix and section 45 of Chandrasekhar 1960), 
i.e. Thomson scattering of unpolarized radiation. This limb 




0.2 0.-1 0.6 O.B 



Figure 4. Angular distribution of total intensity / emerging from 
atmospheres with q = 0.2, q = 0.8, and g = 1; <5 = 0, 2, and 100; 
and all radiation sources at infinite optical depth. In the q = 0.2 
and q = 1 cases the limb darkening is virtually independent of S, 
and the three curves in each case lie nearly on top of each other. In 
the q = 0.8 case, 5 increases from bottom to top. All curves were 
calculated using our Feautrier code. Note that for small q, the 
limb darkening becomes very large, in agreement with equation 

darkening law turns out to be very close to the pure elec- 
tron scattering case with polarization effects. As shown by 
the q — 0.8 curves in figure 4, modest absorption causes 
Faraday rotation to have a more substantial effect on the 
limb darkening, although it is still small (cf. the q = 0.8 
case in figure 4). This is because the polarization is greater 
than for the q = 1 case. The contribution of Q to the inten- 
sity source function is therefore larger, so as the magnetic 
field depolarizes, the intensity source function is modified 
more than for the q — 1 case. For large absorption opacity, 
e.g. the q — 0.2 case shown in figure 4, the effects of Faraday 
rotation are very small. 

Figure 5 shows the depolarizing effects of Faraday rota- 
tion for all values of q. Overall, as q decreases below unity, 
Faraday rotation is at first more effective in depolarizing the 
radiation field than at q — 1. The polarized source function 
has a greater contribution from Q relative to / than it did 
in the q = 1 case. The limb darkening does not change much 
with 5, so as the Faraday rotation is added and the polar- 
ization is reduced (except for /i near 0), the polarized source 
function for fi — decreases more rapidly than in the q = 1 
case. However, as q gets below around 0.4 for the particu- 
lar viewing angle shown in figure 5, Faraday rotation starts 
to become less effective in depolarizing the radiation field 
compared to the 5=1 case because of the diminishing elec- 
tron column density down to unit optical depth. As g — » 0, 
Faraday rotation has zero effect as we noted in our analytic 
solution above. 



© 1996 RAS, MNRAS 000, 



6 E. Agol et al. 





Figure 5. Polarization as a function of q along the /i = 0.452 
line of sight for various values of Faraday rotation parameter <5. 
In (a) we show the actual polarization, while in (b) we show the 
ratio of the polarization to that of the 5 = (unmagnetized) case. 
Prom top to bottom in both figures, the curves represent 5 = 0, 
2, 5, 10, and 100. All the results shown were calculated using our 
Peautrier code. 



3.2 Case 2: Linear Thermal Source Function 

We now consider a distribution of sources in the atmosphere 
such that the thermal source function depends linearly on 
optical depth, 



S(r) = S(0)(l + /3r), 



(22) 



where S{Q) and P are constants. In the q ^ 1 limit, this 
problem reduces again to the pure electron scattering case 
considered in paper I.^ We use the diffusion approximation 
to apply a lower boundary condition at sufficiently high op- 
tical depth Tmax so that the results are independent of Tmax. 
Apart from the uninteresting normalization factor, the radi- 
ation field emerging from this atmosphere now depends on 
three parameters: 5, g, and (i. 

For atmospheres with very small scattering opacities 
(g — > 0), the radiative transfer equation may be solved per- 
turbatively. To lowest order in q, the total intensity from the 
atmosphere is given by the Eddington-Barbier relation, 



/(0,m) = S(0)(1 + /3m) + O(q), 
and the polarization is given by 

3(1 -/i' 



(23) 



16(1+ /3m) 



+o(q' 



(24) 



where 



t This is true provided the thermal source function remains finite. 
By g — > 1, we mean that riecry 3> k and that the scattering source 
function is much larger than the thermal source function. 



Figure 6. Polarization as a function of q along the fj, = 0.452 
line of sight for various values of source function gradient /3 in 
an unmagnetized (5 = 0) atmosphere. Prom bottom to top, the 
curves represent the results of our Peautrier calculations for /3 = 0, 
0.5, 1, 2, 5 , and oo. Points represent the numerical calculations of 
Loskutov & Sobolev (1979), interpolated to ^l = 0.452. The dotted 
lines are the results of the analytic formula for small g, equation 
([2^. All curves approach the Chandrasekhar (1960) value for this 
viewing angle as q ^ 1. 



<1>(m) = f - 3m + (3m' - 1) In 



1 + i 



(25) 



These results are identical to those obtained by Gnedin & 
Silant'ev (1978) for unmagnetized atmospheres in the g ^ 
limit, and this is because the effects of Faraday rotation are 
of order q^ in this limit. Physically, the ratio of scattered to 
thermal (unpolarized) intensity is of order q, which is why 
the polarization is also of this order. (This is in marked con- 
trast to the previous case we considered where the thermal 
sources were all at infinite optical depth. Here the presence 
of a nonzero thermal source function ensures that the polar- 
ization vanishes as q 0.) The amount by which Faraday 
rotation can further depolarize the radiation is also of order 
q, because this is the factor by which the electron column 
density (which does the rotation) down to unit optical depth 
is reduced. Hence we immediately conclude again that Fara- 
day rotation has negligible effect on the polarization (which 
is already small) as g — > 0. 

Figure 6 shows the polarization viewed along /i = 0.452 
for various values of /3 and g for an unmagnetized atmo- 
sphere]^ Also shown are the numerical results of Losku- 
tov & Sobolev (1979), which are again in good agreement 
with ours. Negative values of P in this figure represent cases 
where the polarization plane is perpendicular to the plane of 
the atmosphere. We call this the "Nagirner effect", after the 
person who first noted that absorption opacity can produce 
this (cf. Gnedin & Silant'ev 1978). It is possible that this 



t The /3 = 00 case corresponds to 3{t) oc r. 

© 1996 RAS, MNRAS 000, 



Effects of absorption and Faraday rotation on accretion disk polarization 7 




1=1 



6^ 



-4 -a a 

log(T) 

Figure 7. Polarized source function vs. optical depth along the 
11 = 0.211 line of sight for /3 = and q = 0.8. The solid curve 
is the total source function, while the dashed and dotted curves 
are the contributions from various angles in the four-stream ap- 
proximation (negative fi is downwards). The plot of Sjq(t) for 
the fj, = 0.789 line of sight is very similar. 

efltect could explain the fact that the observed polarization 
in type I AGN is parallel to the radio axis. We find that 
negative polarization is present for some fi and q if and only 
if /3 < (6 - 81n2)/(81n2 - 5) = 0.834. This is consistent 
with Gnedin & Silant'ev's (1978) g — » result, because as (3 
drops below the critical value, the polarization first becomes 
negative for small q. 

The Nagirner effect therefore arises for sufficiently fiat 
thermal source functions. It is useful, however, to examine 
its origin a little more closely by considering the depth de- 
pendence of the total source function. If there is no magnetic 
field, then the outgoing radiation (0 < p < 1) can be for- 
mally expressed from equation (^^ as 

1(0, m)= / 'S{t,fi)e~'^^-. (26) 
Jo 

Hence if £sQ(r, /i) < over some range of optical depths, 
then the outgoing radiation can be negatively polarized for 
that particular viewing angle. The solid curve in figure 7 
shows the depth dependence of Sjq(t") for l3 — 0, q = 0.8, 
and fi = 0.211. Around t — 1, 3q is negative, but for low 
r it becomes positive. The reason for this can be seen by 
looking at the contribution to the polarized source function 
from radiation coming from different directions in an atmo- 
sphere with constant source function. The broken curves in 
figure 7 show the contribution to the source function from 
different angles in a four-stream calculation of the radiation 
field, i.e. where the radiation is calculated at four angles for 
the purposes of computing quadratures. Near r = 0, there 
is no downgoing radiation {fi < 0), so the only contribu- 
tion is from upgoing radiation. The limb darkening causes 
the source function polarization to be positive, since the 
near-vertical radiation (which has net positive polarization 
when scattered) is stronger than the near-horizontal radia- 




Figure 8. Cartoon showing the reason why the polarized source 
function switches from positive to negative to zero with increas- 
ing T. The length of the radial arrows represents the strength of 
radiation coming from different angles. The arrows in the centre 
represent the strength of the negative and positive contributions 
to the polarized scattering source function (negative is vertical, 
positive is horizontal). The sum of the polarizations is indicated 
on the right, along with the optical depth of each layer. 



tion (which has net negative polarization when scattered). 
Near r = 1, the limb darkening is weaker. Here, there be- 
gins to be a significant contribution to the source function 
from downward radiation which is produced in the layers 
above r = 1. The near-horizontal radiation is stronger than 
the near-vertical since there is more atmosphere emitting 
from smaller \^\. This leads to a net negative polarization 
at T = 1. Figures 8 and 9 illustrate this effect. (Note that the 
polarization is low, so most of the contribution to comes 
from /. Thus, the corresponding difference in the contribu- 
tion to 9q from different angles is due mostly to the differ- 
ence in the intensities from different angles.) When the ther- 
mal source function has a steep vertical gradient, the limb 
darkening is stronger. The contribution from the downward 
going radiation is less than that of the upward radiation, 
and the polarization source function is always positive. This 
is why flat thermal source function gradients are required 
for the Nagirner effect to be present. 

Figure 10 shows the polarization as a function of view- 
ing angle for l3 = 1 and (a) q = 0.8 and (b) q = 0.2. This 
figure should be compared with figure 2 above. Faraday ro- 
tation has a much smaller effect on the polarization at /i = 
when a nonzero thermal source function is present. This is 
because the polarization is low and limb darkening is more 
important in determining the polarized source function. In 
addition, the presence of the thermal source function im- 
plies that the overall polarization vanishes as </ — > 0. Fara- 
day rotation has even less effect on the polarization, which 
is already small, as g ^ 0. For example, in the q = 0.2 case 
shown in figure 10(b), the S = 0, 2, and 5 curves overlap 
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Figure 9. At wavelengths where the Nagirner effect is present, 
the polarization vector lies in the plane of the line of sight and 
the vertical at r = 1, while it is perpendicular near r = 0. Thus, 
the polarization ends up being negative at some viewing angles. 
However, when a magnetic field is added, Faraday rotation causes 
the polarization near t = 1 to be rotated and depolarized, while 
the radiation near r = is rotated very little, so the polarization 
ends up being nearly perpendicular. 

because the effects of Faraday rotation are reduced by an 
additional factor of q as discussed above. 

Figure 11 shows the polarization in the (3 = 1 and oo 
cases as a function of q for various 5. Notice again that as q 
becomes small, the effect of the Faraday rotation decreases, 
and all curves approach the 5 = case. Note however that 
this approach is much smoother than in the previous case 
where all the sources were at infinite optical depth (cf. fig- 
ure 5). For moderate absorption opacity {q slightly less than 
unity), the depolarizing effects of Faraday rotation are again 
enhanced over the pure scattering problem for the P = oo 
case shown. This is the same effect as in section 3.1, and 
again arises because the absorption opacity on its own in- 
creases the polarization. If the thermal source function gra- 
dient is not so steep so that modest absorption opacity de- 
creases the polarization, then the effects of Faraday rotation 
are reduced as q drops below unity (cf. the P = 1 case in 
figure 11). 

So far we have discussed Faraday rotation as an agent 
for depolarizing the radiation field. However, when 13 is small 
and the Nagirner effect is present, it is possible for the mag- 
netic field to increase the polarization. Figure 12 shows the 
normalized Stokes parameters for g = 0.9 and /3 = 0.25 for 
various 5, and Q/I with the 5 = case subtracted. With 
no magnetic field, the Sjq's from different depths partially 
cancel, making the polarization negative for > 0.38, but 
positive for smaller /i. 

Since = 0, the only contribution to the polarization 
is from 9q, which for q = 0.9 and (3 = 0.25 is similar to the 





(a) : 
- 




(b) : 







0.2 0.-1 6 0.8 1 

M 

Figure 10. Polarization as a function of viewing angle for (a) 
q = 0.8 and (b) g = 0.2 atmospheres with various values of 5 
and a linear source function with (3 = 1. From top to bottom, 
the curves represent the results of our Feautrier calculations for 
5 = 0,2, 5, 50, and 100. 

source function plotted in figure 7. Because the polarization 
is low, 9q does not change much with 5 here. For this case, 
3q < for r > 0.25 at all /i. Thus, the positive polarized 
source function comes from a small range of electron scat- 
tering depths, so the Faraday depolarization is insignificant 
for (5 < 4. However, the negative polarized source function 
comes from a larger range of optical depths, causing Fara- 
day depolarization for 5 ~ 1, reducing the magnitude of the 
negative polarization so that the outgoing polarized fiux is 
increased. ^ For very large 5, even the radiation from small r 
will be Faraday rotated significantly, so the polarization will 
still be approximately horizontal, but will eventually start 
to decrease in magnitude again, as seen in figure 12. The 
effect of increasing the polarization with increased magnetic 
field in a semi-infinite atmosphere is only present when there 
is absorption opacity present and a shallow source function 
gradient. 



4 POLARIZATION FROM REALISTIC 

ATMOSPHERES IN AN ACCRETION DISK 

In a previous paper (Blaes & Agol 1996), we calculated the 
structure of local, static, plane-parallel atmospheres using 

§ Strictly speaking, "positive" and "negative" polarization define 
the orientation of the plane of polarization only when there is no 
magnetic field present, because there are only two possible orien- 
tations. The presence of the magnetic field breaks the azimuthal 
symmetry and allows the plane of polarization to be at an arbi- 
trary angle with respect to the vertical/line of sight plane. We 
use these terms here in reference to the unmagnetized case to 
show how Faraday rotation acts on the different orientations of 
the polarization at different depths in the atmosphere. 
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Figure 11. The effect of Faraday rotation on the /3= 1 (dashed) 
and oo (sohd) cases for the /i = 0.452 hne of sight. In (a) we 
show the actual polarization, while in (b) we show the ratio of 
the polarization to that of the 5 = (unmagnetized) case. From 
top to bottom in both figures, the curves represent the results of 
our Feautrier calculations for 5 = 0, 2, and 5. 
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Figure 12. The top panel shows P and Q/I, and the bottom 
panel — Q{0))/I and U/I as 5 increases in an atmosphere 

with q = 0.9 and f3 = 0.25. Note that at large 5 and small fi, Q/I 
decreases relative to the zero magnetic field case, which is because 
the Faraday rotation is significant for small optical depths when 
S is large enough. 



Figure 13. Stokes parameters in an atmosphere including Fara- 
day rotation for T^ff = 20,000 K and g = 190 cm s~^ includ- 
ing non-LTE effects for two hydrogen levels. The curves depicted 
are for a, /i = 0.452 line of sight, with various magnetic field 
strengths. The top panel shows the outgoing total intensity spec- 
trum Ix = vl\Iv The bottom two panels show the percent po- 
larization for the Q// and Vjl Stokes parameters. Qjl peaks 
at 10 per cent and {/// peaks at -4 per cent at about 500A for 
B = 750 G. 

the complete linearization technique, neglecting the effects 
of any magnetic field .|^ Here we use some of these atmo- 
spheres to calculate the radiative transfer in the presence 
of a constant vertical magnetic field in the atmosphere ac- 
cording to equation ([T4[). Our atmosphere solutions include 
hydrogen bound-free and free-free opacities as well as elec- 
tron scattering opacity. Non-LTE effects in the n = 1 and 2 
levels of hydrogen are included. The results with magnetic 
fields of different strengths in atmospheres with {Teff, g) of 
(2 X 10'' K, 190 cm s'^), (4.5 x lO" K, 4 x 10^ cm s^^), 
and (10^ K, 9.5 x lO"* cm s"^) are plotted in fi gures 13-15, 
respectively. We also show the variation of g at r = 1 with 
wavelength in figure 16 for these three atmospheres. In all 
cases the magnetic field does not significantly affect the total 
intensity spectrum I\. 

As shown in figure 16, q is very small just blueward 
of the Lyman edge, so the Faraday rotation does not affect 
the polarization very much in this region of the spectrum 
unless 5 is very large. This is illustrated in figures 13-15. 
In the case shown in figure 13, the source function is shal- 
low just blueward of the Balmer edge. Hence 9q < near 
r = 1 at these wavelengths. Here the increasing magnetic 
field causes an increase in polarization, as described in sec- 
tion 3. Figure 13 also shows a dramatic reduction in the 
difference in polarization across the Lyman edge. This is 
due to the fact that redward of the edge, is large, so 



^ In particular, we ignore the contribution of magnetic field pres- 
sure on hydrostatic equilibrium, which may be quite important 
for equipartion fields. 



© 1996 RAS, MNRAS 000, 0-|l| 



10 E. Agol et al. 



D 





Figure 14. Polarization in an atmosphere including Faraday ro- 
tation for T^ff = 45,000K, g = 4000 cm s'^, and ^l = 0.55 and 
various magnetic field strengths. 



MA) 

Figure 15. Polarization near the Lyman edge in an atmo- 
sphere including Faraday rotation for T^f f = 100, 000 K, g = 
95, 000 cm s"'^, and ^ = 0.55 and various magnetic field strengths. 



the Faraday depolarization is strong. Blueward of the edge, 
however, is small, so the Faraday depolarization is weak. 
Thus, as the magnetic field increases, the polarization de- 
creases faster red ward than blueward of the edge. In some 
cases the polarization blueward of the edge is larger than 
redward of the edge, as shown in figure 15 for B — 2500 G. 

The Balmer edge in figure 14 has a negative polariza- 
tion, which becomes more negative redward of the edge for 
B = 0. However, as the magnetic field increases, the polar- 
ization first becomes positive, and increases from blue to red 
across the edge. For larger B, the polarization is reduced, de- 
creasing from blue to red. This behavior can be understood 
in terms of the physics discussed in section 3. The thermal 
source functions are very flat on both sides of the Balmer 
edge, but q is larger redward of the edge (q ~ 0.8 at r = 1) 
than blueward of the edge {q ~ 0.57 at r = 1). Hence the 
Nagirner effect is stronger (cf. figure 6). Since Faraday depo- 
larization is stronger for larger q, the polarization increases 
to the red across the edge when there is a 20 G magnetic field 
(5 = 2.1 at 3648A). However, when S is larger, the depolar- 
ization is more rapid for larger q, so then the polarization 
decreases to the red across the edge when B — 100 G, or 
S = 10.5 at 3648A. 

If the magnetic fleld were randomly oriented, the polar- 
ization feature may also be decreased at larger inclination 
angles for a large magnetic fleld. 

The results of this section are meant to give a flavor 
of the complexities that may result in the polarized spec- 
trum emerging from a realistic atmosphere. The full accre- 
tion disk spectrum requires an integration over many such 
atmospheres in different physical conditions representing the 
disk at different radii. 
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Figure 16. The ratio of scattering opacity to total opacity at 
T = 1 for the three atmospheres displayed in figures 13-15. The 
solid, dotted, and dashed curves represent qx for T^ff= 20,000K, 
45,000K, and 100,000K atmospheres respectively. 

5 CONCLUSIONS 

Using a combination of numerical calculations and analytic 
arguments, we have extended our previous study of Faraday 
rotation in accretion disk atmospheres to include the interac- 
tion with absorption opacity. Along the way we have clarified 
the role that each of these two effects play separately. Fara- 
day rotation in a pure electron scattering atmosphere acts 
to depolarize the radiation field by rotating the polarization 
vectors of photons scattered from different depths. 
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Absorption opacity in an unmagnetized atmosphere can 
increase or decrease the polarization depending on the be- 
havior of the thermal source function. If the thermal source 
function is zero except at great depth, absorption opacity 
alone always increases the polarization. In the more usual 
case of nonvanishing thermal source function, the effect of 
absorption opacity depends on the source function gradient. 
If the thermal source function increases steeply with depth, 
absorption opacity can increase the polarization over the 
pure electron scattering case. On the other hand if the source 
function increases slowly with depth, or even decreases, then 
absorption opacity can flip the plane of polarization to be 
in the vertical/line of sight plane. While these results were 
known from previous work (Gnedin & Silant'ev 1978, Losku- 
tov & Sobolev 1979), we have presented a novel physical 
interpretation in terms of the behavior of the total source 
function Ss. Quite generally this source function locally pro- 
duces polarization which is parallel to the atmosphere plane 
at low optical depths. At optical depths around unity, how- 
ever, the polarization can be parallel or perpendicular to the 
atmosphere plane depending on the thermal source function 
gradient. 

When Faraday rotation is combined with absorption 
opacity in a scattering atmosphere, a number of effects can 
occur. First, if absorption dominates scattering, then the 
electron column density along a photon mean free path is 
small and Faraday rotation only has a small effect on the po- 
larization of the radiation field. On the other hand, if mod- 
est absorption opacity alone increases the polarization, then 
the depolarizing eff'ects of Faraday rotation are enhanced 
compared to the pure scattering case, at least for a verti- 
cal magnetic field. On the whole Faraday rotation generally 
acts to depolarize the radiation field, as in the pure scat- 
tering case. However, when the Nagirner effect is present at 
certain photon wavelengths, then Faraday rotation can actu- 
ally increase the emerging polarization by depolarizing the 
deeper radiation field which has a perpendicular orientation 
to the radiation which is scattered from shallower depths. 

While we have shown that these effects can all occur in 
simple toy model atmospheres, we have also demonstrated 
that they are present in more realistic atmospheres. We have 
also shown how the polarized radiative transfer can be com- 
puted in a straightforward manner by a simple extension 
of the Feautrier method to incorporate a vertical magnetic 
field. This numerical method can be used to integrate the 
radiation field produced at different annuli in the accretion 
disk to calculate the total observed radiation field. As we 
noted in paper I, the largest uncertainty in applying such 
calculations to the observed data is the variation of mag- 
netic field strength with disk radius, and its covering fac- 
tor on the disk photosphere. (It might be possible to get 
a handle on the latter by combining ultraviolet and X-ray 
observations to determine the "patchiness" of the corona; 
cf. Haardt, Maraschi, & Ghisellini 1994. This assumes that 
such patches, if real, are magnetized active regions similar 
to those on the sun.) The field topology will of course also 
not be vertical, but following our Monte Carlo work of paper 
I, we expect magnetic fields of random orientation to have 
qualitatively (and perhaps quantitatively to some extent) 
similar effects. 
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APPENDIX A: 

We reproduce here Silant'ev's (1979) analytic calculation 
of the polarization of radiation emerging from an optically 
thick, pure electron scattering atmosphere in the S ^ oo 
limit. 

Silant'ev's calculation uses a radiation density matrix 
formalism, which is applicable to magnetic fields with arbi- 
trary orientation. For a vertical magnetic field, however, it 
is somewhat simpler to use the radiative transfer equation 
as we have formulated it in section 2. We therefore proceed 
to show how Silant'ev's result can be derived within this 
formalism. 

Define a complex intensity column vector 
/ 7/21/2 . 
I'= i~Q + iU)/2 . 
\{-Q-iU)/2) 



(Al) 

Then, for q — 1, the equation of transfer ( |l^ may be written 

riV 1 fi 

M|^=Ml'-iy VP'(m,m')I'(m'), (A2) 

where 

/I \ 

M = 1 =F j'JM , (A3) 

\0 l±i5nJ 

P'ij{li,li') = ai{ii)aj{fi') +5iiSji, (A4) 

and 

a(M) = (-P2(p)/2l/^ 3(1 - m')/4, 3(1 - m')/4) . (A5) 
As in section 2.2, the upper and lower signs in equation 



(A5) refer to upward and downward magnetic field direc- 
tions, respectively. We have again suppressed the frequency 
subscript because there is no frequency redistribution in this 
problem. 
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Using a principle of invariance for vectors (of. chapter 4 
in Chandrasekhar 1960), the radiation field emerging from 
the atmosphere (0 < < 1) can be expressed in terms of a 
scattering matrix S: 

/I ON 



\0 



P'(^i,^')I'(0,m')V 



M MO / 

= P'(/.,Mo) + i /'p'(M,M')S(/i',Mo)^ + 
^ Jo ^ 

i/ S(/i,/i')P'(M',Mo)^ + 
^ JO ^ 

S(/i, /i')P'(^i', ^")S(/i", Mo)^ 



(A6) 



+ i / / S(m,m')P'(m',m")I'(0,m")^ 
* Jo Jo 

The scattering matrix satisfies the following integral equa- 
tion: 




(A7) 



(A8) 



Note that the symmetry of equations (A4) and (A7) implies 
that S(//, ^o) = S(^o,m)- 
Define 

/ 1 l±iSC l^iS( \ 
7(/i, fio) = ^l + ^10 i I T iSC, 1 It 2i<SC , (A9) 
Vl±i(5C l±2i5C 1 / 

(where = fj,fj,o/{^i + /io)) and a matrix T(/i, /to) by letting 

/i/iO 



SiJ{^^,^^o) 



-Tij{fi,fio) 



(AlO) 



7i:i(M,Mo)" 

where no summation over the indices i and j is implied. In 
the limit of large 5, l/"fij{n, no) = Sij/{fi + fio), provided 
both fj, and /to are nonzero. If we make the ansatz that all 
elements of the matrix T remain finite for all values of /i 
and fio, then in this limit equation (A7) simplifies to 



TijifJ., fJ-o) = Hi{fi)Hj{fj.o) + K,{p)Kj{fio)- 
Here 

H^{^i) = ai{ii) + - / ^ , d^l , 



Ki{p) = Sii 



1 + ^ f ^M^d/.' 
/i + /i' 



(All) 



(A12) 



(A13) 



and we are considering both /i and /to to be nonzero. Since 
12 (/t) = 03 (m)) the equations for T22 and T33 are identical. 
Let i/2(/i) = Hsifi) = 1(1 - m^)<^(m)- Then (/>(/t) satisfies 
the nonlinear equation 



M = 1 + 



32 



to M + M' 
The equation for Tn may be written 



d/t'. 



Til (/t, /to) 



1 8 ~ ~ 

-V'(/i)V'(/'o) + -^'i'M<f>{^J■o) 



(A14) 



(A15) 



where 



i/;(/i) = 3 - /i^ + 



/i Iii(m^mO(3_V) 
/i + /i' 



d/t' 



and 



^Jo ^ + 



(A16) 



(A17) 



Apart from multiplicative factors, equation (A15) is identi- 
cal to the equation for the first term of the scattering func- 
tion S^'^^fi, fio) in the problem of diffuse refiection from a 
scattering medium described by Rayleigh's phase function 
(cf. equation 6 in section 44 of Chandrasekhar 1960). 

The scattering matrix elements in the limit of large S 
are therefore given by 

S'ii(/i,/to) = I jj-(/t)jj-(/to)[3-c(/t + /to)+/^MojlA18) 

5 /t + /to 



5*22 (a*, Mo) 



= S'33(/i,/io) 
^ 9 flflQ 

16 /t + /to 



-Atg)0(/i)0(/io) (A19) 



where //(/t) satisfies the integral equation. 



= 1 + j^fiHifi) 



-H{n')dfi', 



and 



Jo H{fj,)fj,dfi 



(A20) 



(A21) 



The off-diagonal components of S(/t, /to) are negligible (~ 
5~^). Although these expressions were derived assuming 
both /i and no are nonzero, they remain valid even if this 
is not the case. This is because they then vanish provided 
T(/i,/io) is finite, which can easily be shown to be true. 

Given this solution for S(/t, /io), we can calculate the 
outgoing radiation field from equation (A6). We immedi- 
ately conclude from the matrix multiplier in this equation 
that Q and U vanish as 5 ^ 00 for nonzero /i. The intensity 
is given by 



AO,At) = i 



P'ii(/i,/i')/(0,/i')d/i' 
d/t' 



Sn(/i,/.')P'n(/i',/i")/(0,/.")^d/t", (A22) 



which has solution 

FH{n) 



7(0, /i) = 



27T fi'H{fi')dfi'' 



(A23) 



where F is the fiux emerging from the atmosphere. This is 
the same limb darkening law as that of radiation emerg- 
ing from a scattering medium described by Rayleigh's phase 
function (cf. section 45 of Chandrasekhar 1960), i.e. an elec- 
tron scattering medium in which the radiation is everywhere 
treated as unpolarized. This makes physical sense in the 
limit of infinite Faraday rotation. 

Since the Q' and U' Stokes parameters are of minimum 
order (for fi not equal to 0), we can ignore the contri- 
bution of Q' and U' in the integrals on the right hand side 
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of equation (|A^). This then gives: 

»i 



Q' 



1 



1 — iSfi 



1 

+ 4 
and 
U' 



dfi' 



522(^1, n')Pii{n', f^")l'iO, ^i") — dfj," 



(A24) 



(A25) 



Performing these integrals, using equation (42), and solving 
for Q and U, we get: 



Q 



3(l-/i^MM) 
8(1 + (5V^) 



and 



U 



3S^l{l - ^^)0(m) 
8(1 + <5V^) 



(A26) 



(A27) 



which are the same as the results obtained by Silant'ev 
(1979). The polarization at = is independent of 5, even 
for 5 = oo, when the polarization is zero everywhere except 
H — where it is equal to 9.137 per cent. We completely re- 
tract our claim in paper I that Silant'ev's treatment breaks 
down near ^ = 0. This was based on inaccuracies in our 
Monte Carlo results at that time. 
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